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ABSTRACT 



One of the principal scientific objectives of the upcoming Constellation- X 
mission is to attempt to map the inner regions of accretion disks around black 
holes in Seyfert galaxies by reverberation mapping of the Fe Ka fluorescence 
line. This area of the disk is likely radiation pressure dominated and subject to 
various dynamical instabilities. Here, we show that density inhomogeneities in 
the disk atmosphere resulting from the photon bubble instability (PBI) can cause 
rapid changes in the X-ray reflection features, even when the illuminating flux 
is constant. Using a simulation of the development of the PBI, we find that, for 
the disk parameters chosen, the Fe Ka and O VIII Lya lines vary on timescales 
as short as a few hundredths of an orbital time. In response to the changes 
in accretion disk structure, the Fe Ka equivalent width (EW) shows variations 
as large as ~100 eV. The magnitude and direction (positive or negative) of the 
changes depends on the ionization state of the atmosphere. The largest changes 
are found when the disk is moderately ionized. The O VIII EW varies by tens 
of eV, as well as exhibiting plenty of rapid, low-amplitude changes. This effect 
provides a natural explanation for some observed instances of short timescale 
Fe Ka variability which was uncorrelated with the continuum (e.g., Mrk 841). 
New predictions for Fe Ka reverberation mapping should be made which include 
the effects of this accretion disk driven line variability and a variable ionization 
state. Reflection spectra averaged over the evolution of the instability are well 
fit by constant density models in the 2-10 keV region. 
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1. Introduction 

Shortly after the prediction of relativistically broadened iron Ka fluorescence lines in the 
X-ray spectra of accreting black holes (Fabian et al. 1989), it was realized that if the lines were 
variable then constraints on the emission geometry and basic black hole parameters could 
be determined (Stella 1990; Matt & Perola 1992). These ideas matured into the prospect of 
iron line 'reverberation mapping', which, given a series of short time- resolved X-ray spectra 
following a significant flare in the continuum, could accurately derive such fundamental 
parameters as the black hole mass, spin, and location of the X-ray source (s) (Reynolds et 
al. 1999; Young & Reynolds 2000; Ruszkowski 2000; Goyder & Lasenby 2004). The promise 
of this technique is such that iron line reverberation mapping of nearby Seyfert Is is one of 
the principal science goals of the upcoming Constellation- X observatory. 

In recent years it has become obvious that reverberation mapping will be more chal- 
lenging than perhaps was previously thought. The XMM-Newton telescope has found only 
a handful of relativistically broadened Fe Ka lines, which indicates that these components 
may be in general quite weak and low contrast features. Moreover, the best example of a rel- 
ativistic line (observed from MCG-6-30-15; Fabian et al. 2002) exhibits less variability than 
the X-ray continuum and shows no correlation with the observed broadband flux (Vaughan 
& Edelson 2001; Vaughan & Fabian 2004, but see Ponti et al. 2004). These properties seem 
to apply on long timescales in a wide sample of bright Seyfert 1 galaxies (Weaver, Gelbord 
& Yaqoob 2001; Markowitz, Edelson & Vaughan 2003), but there also have been instances of 
rapid iron line variation with only a slight continuum change (Mrk 841; Petrucci et al. 2002), 
and where the line at least partially follows the continuum (e.g., Nandra et al. 1997; Iwasawa 
et al. 2004). Clearly, there is a wide range of possible Fe Ka variability characteristics. 

The problem these observations cause with the standard reverberation picture is with 
the underlying assumption that variations in the Fe Ka line are directly related to changes 
in the observed continuum. The observed behavior shows that this assumption cannot 
be valid at all times, and that there are other effects which are equally important to the 
production of the Fe Ka line. One possible perturbation to the previous reverberation 
calculations is photoionization of the accretion disk surface by the incident X-rays (the work 
by Reynolds et al. 1999, e.g., assumed only neutral reflection from the disk outside of the 
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innermost stable circular orbit). If the disk surface is ionized to a certain level, then the 
Fe Ka equivalent width (EW) becomes anti-correlated with the continuum, in contrast to 
the positive correlation expected if the disk remains neutral (Nayakshin 2000; Ballantyne, 
Fabian & Ross 2002; Ballantyne & Ross 2002). It is also possible that the X-ray continuum 
incident on the disk is different from the observed one. For example, light-bending effects 
could enhance the illumination (Miniutti & Fabian 2004), while an outflowing corona could 
diminish it (Beloborodov 1999). While these effects may all be operating at some level, they 
are still connected with the X-ray continuum, and therefore a relation (positive or negative; 
strong or weak) between the Fe Ka line flux and the incident continuum is expected. In 
this paper, we outline a process by which the Fe Ka line may vary independently of the 
illuminating X-ray flux. 

For systems accreting at greater than 0.3% of the Eddington rate, the innermost annuli 
of the accretion disk, where the relativistic portion of the iron Ka line is produced, are likely 
radiation pressure dominated (Shakura & Sunyaev 1973). Numerical results indicate that 
turbulence driven by the magneto-rotational instability (MRI; Balbus & Hawley 1991) may 
involve large density fluctuations in radiation-dominated disks if magnetic pressure exceeds 
gas pressure and photons diffuse from compressed regions (Turner et al. 2003; Turner 2004). 
Large density variations may also be caused by the photon bubble instability (PBI; Arons 
1992; Gammie 1998), which may develop into a series of shocks propagating through the 
plasma (Begelman 2001). These results suggest that radiation-dominated disks have time- 
varying density inhomogeneities on scales smaller than the disk thickness. In the first paper 
of this series (Ballantyne, Turner & Blaes 2004, hereafter Paper I), we showed that density 
variations in the outermost few Thomson depths can lead to reflection spectra differing 
from those of uniform slabs. Here, we continue our investigation of X-ray reflection from 
inhomogeneous accretion disks by computing the evolution of the reflection spectrum as 
the disk structure is changing due to the PBI. As is shown below, the changing density 
structure in the disk can affect the EW and line flux of the predicted emission lines, despite 
the constant illuminating continuum. 

The next section details the setup and conditions of the PBI model and reflection calcu- 
lations. We then present the results of the reflection models, both as a function of time and 
time-averaged (§ 3). The results are then discussed in § 4 before we summarize and draw 
our conclusions in § 5. 
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2. Models: Setup and Calculations 

An important timescale in X-ray reflection is the time for the accretion disk to re- 
establish hydrostatic balance following a change in illumination (Nayakshin & Kazanas 2002; 
Collin et al. 2003). The adjustment time is about equal to the orbital period (Frank, King & 
Raine 2002). Turbulence inside the disk results from the MRI and also has a characteristic 
overturn time similar to the orbital period (Balbus & Hawley 1998). In radiation-dominated 
disks, the PBI may grow faster than the orbital frequency (Blaes & Socrates 2003), and 
could lead to structures that vary over a small fraction of an orbit. Below we examine time 
variations in the X-rays reflected from structures produced by the growth of photon bubbles. 
The variations are computed in two steps, as in Paper I. First, density distributions resulting 
from the instability are found using a two-dimensional, time-dependent numerical radiation 
magnetohydrodynamical (MHD) calculation. The density profile along a single ray passing 
vertically through the disk is used as the input for the second step. We assume X-rays from 
a diffuse external source strike the disk, and calculate the reflected spectrum by solving the 
coupled equations of radiative transfer and ionization balance in one spatial dimension. The 
effects of changes in the density profile are tracked by calculating the spectrum 100 times 
per orbital period. 

The calculation of the growth of photon bubbles is similar to that described in Paper I. A 
narrow annulus of accretion disk in orbit around a black hole of M BH = 1O 8 M at a distance 
Rs = 20 Schwarzschild radii (= A0GM B u/c 2 ) is modeled in locally Cartesian coordinates in 
two dimensions, with symmetry assumed along the orbit. The frequency-averaged equations 
of radiation MHD are integrated using the ZEUS code (Stone & Norman 1992a,b) with 
its flux-limited radiation diffusion module (Turner & Stone 2001), on a grid of 128 x 768 
uniformly-spaced zones. The initial condition is a Shakura-Sunyaev model with accretion 
rate 10% of the Eddington limit for 10% radiative efficiency, and accretion stress parameter 
a = 0.06. The total Thomson optical depth is 9200. To this is added a uniform magnetic field 
with pressure 1% of the midplane radiation pressure, inclined 12° from horizontal. The initial 
state is disturbed slightly by applying random density perturbations of up to 1% in each grid 
zone. The domain extends 1.2 Shakura-Sunyaev semi-thicknesses either side of the midplane 
and has width one-sixth the height. The horizontal boundaries are periodic, and the vertical 
boundaries allow gas, radiation and magnetic fields to flow out but not in. Differential 
rotation is neglected, so there is no magneto-rotational instability and no external energy 
source. The radiation initially present escapes through the top and bottom boundaries, and 
the gas cools over time. The flux of radiation energy leaving each face of the disk initially is 
-Fdisk = 7.7 x 10 13 erg cm~ 2 s -1 . During the first 0.8 orbits of the calculation, photon bubbles 
grow to non-linear amplitudes and develop into trains of shocks. Photons then diffuse faster 
through the low-density regions between the shocks and the disk flux is several times greater 
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than initially. By 2.35 orbits, three-quarters of the initial radiation energy is lost, and the 
horizontally- averaged disk thickness approximately halves. A density floor of 1% of the initial 
midplane value or 9.93 x 10~ 12 g cm" 3 is imposed during the calculation so that results are 
obtained with a reasonable amount of computer time. The floor is applied in the regions 
between shocks, reducing the overall density contrast. The dynamical time for this radius is 
Wn ~ (R 3 /GM BR y/ 2 ~ V8(R/R s ) 3 / 2 GM B}i /c 3 ~ 1.2 x 10 5 ( J R/20 i? s ) 3/2 (M BH /10 8 M ) s. 
The orbital period is 7.8 x 10 5 seconds, or about nine days. The density profile along a 
ray passing vertically through domain center is stored every 0.01 orbits for use in the X-ray 
reflection calculations. 

Each density profile was trimmed so that the total Thomson depth through the gas 
was about 10. The number of zones in the density cuts varied from a minimum of 23 
to a maximum of 71. Reflection spectra were then computed for each profile using the 
code described by Ross & Fabian (1993) (see also Ballantyne, Ross & Fabian 2001). The 
illuminating continuum was a power-law defined between 1 eV and 100 keV with a photon 
index T = 2. The gas was allowed to come into thermal and ionization balance before the 
reflection spectrum was computed (these timescales are generally much shorter than the 
hydrostatic time; e.g., Nayakshin & Kazanas 2002). The following ions are included in the 
calculations: C V - VII, N VI - VIII, O V - IX, Mg IX - XIII, Si XI - XV and Fe XVI 
- XXVII. Four different illuminating fluxes F x were used for each density profile in order to 
investigate differences due to ionization effects. These fluxes corresponded to Fx/ F disk =0.5, 
1, 4, and 8. Here, Fdi S k corresponds to the initial disk flux, and Fx does not change with 
time. All 944 ionized reflection models converged with no problems. 

Of course, this procedure is not self-consistent, as the heating due to the illuminating 
X-rays is not included in the gas dynamics. Similarly, the internal radiation field in the MHD 
calculation does not change the ionization state of the gas, nor is it included in the outgoing 
spectrum. Unfortunately, time-dependent 2-D and 3-D calculations of an X-ray illuminated 
accretion disk with ionization physics, frequency- dependent radiative transfer and MHD are 
still years away. However, we expect that the omissions described above do not greatly affect 
our results. The internal radiation field has temperatures ~ 10 5 K, so neglecting its ionizing 
power does not affect the high energy spectrum that is the focus of this paper. Ignoring the 
effects of X-ray heating on the gas structure is more problematic, as an optically thick skin 
may form on the disk surface in hydrostatic models with Fx/F disk > 1 (e.g., Nayakshin, 
Kazanas & Kallman 2000). This skin would then scatter and smear out the reflection 
spectrum. However, the return to hydrostatic balance after a pressure perturbation takes 
about an orbit, so the MRI and PBI may cause permanent disequilibrium. The turbulence 
driven by the MRI can lead to density changes in the gas on the orbital period (Turner et 
al. 2003) and can mix irradiated gas in to large optical depths. The fastest photon bubble 
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modes have growth time shorter than the orbital period by the square root of the ratio of 
radiation to gas pressure. This estimate follows from Table 2 of Blaes & Socrates (2003) 
if radiation and magnetic pressures exceed the gas pressure, as in the surface layers in a 
3-D radiation-MHD calculation of a small patch of disk that includes the generation and 
buoyancy of the magnetic fields (Turner 2004). Furthermore in the numerical photon bubble 
calculation described at the start of this section, the instability leads to a radiative flux with 
horizontal component comparable to the vertical. The horizontal flux varies across the face 
of the disk so the atmosphere feels substantial differential stresses. As is seen below, the 
densities change over periods as short as a few percent of an orbit, and the gas is out of 
hydrostatic balance from the time the instability reaches non-linear amplitudes until the end 
of the calculation. These results suggest that like the disk itself, any X-ray heated skin is 
inhomogeneous and time- variable. The effects of irradiation on a turbulent atmosphere may 
be worth investigating with future calculations. 



3. Reflection Results 

3.1. Time-dependent Variations 

To quantify any changes made to the reflection spectra by the evolution of the PBI, we 
calculated the Fe Ka and O VIII Lja EWs from each spectrum. Since the T = 2 power-law 
will also be observed along with the reflection component, the two spectra were summed 
before the EW calculations were performed. The results are shown in Figure 1, with the 
different lines denoting the different values of F x / F disk . The Fe Ka EW spans the range 
from 20-60 eV when Fx/F disk =1 and the line is suppressed due to Auger destruction (Ross, 
Fabian & Brandt 1996), up to ~ 320 eV when Fx/-Fdi S k =8 and a strong ionized line at 
6.7 keV is present. The plots show that there was little variation in the EWs of the emission 
lines over the first orbit, but large changes occurred afterward as the PBI grew in magnitude. 
A very significant variation in the Fe Ka EW occurred between orbits 0.9 and 1.2, where the 
density at the outermost layers of the atmosphere increased by ~ 50%. Despite the small 
increase in density, it greatly altered the Fe Ka EW, particularly when the line originated 
from ionized iron (i.e., when -F x /F disk =4). In that case, the Fe Ka EW dropped by a 
factor of ~ 2.5 due to a decrease of the effective ionization parameter and an increase of the 
continuum absorption. On the other hand, when a neutral 6.4 keV line dominated or when 
the atmosphere was very highly illuminated, the small increase in density between orbits 1.0 
and 1.1 caused much smaller changes in the Fe Ka EW (20-40 eV). 

As mentioned in § 2, a density floor equivalent to 4.25 x 10 12 cm -3 was used in the PBI 
simulation. Tests with a floor artificially lowered by a factor of 10 resulted in a more ionized 
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reflection spectrum, as if F x had been increased and the density was left unchanged. We 
conclude that the density floor will not qualitatively affect the results presented in Fig. f . 

As an illustration of the rapidity of the possible changes to the reflection spectrum, 
Figure 2 plots the model spectra between 0.01 keV and 20 keV for every 0.01 of an orbit 
between orbits 1.6 and 1.69. The illuminating flux is -Fx/-F disk =4, and, as seen in Fig. 1, 
the Fe Ka EW drops and recovers by ~ 60 eV over this range. The inset in each panel 
shows the hydrogen number density nn profile used in the reflection calculation and the 
equilibrium gas temperature found from the model. These plots illustrate that the density 
profile in the outer 10 Thomson depths of the accretion disk can change significantly over as 
little as 0.01 of an orbit. Some of the changes affect the reflection spectra, and some of them 
do not (depending on how strongly the atmosphere is photoionized) . However, in this case, 
a rapid increase in density at the surface of the atmosphere caused a significant variation in 
the strength of the predicted Fe Ka line. 

The second region of large amplitude line variability takes place about a full orbit 
after the first instance, between orbits 1.9 and 2.2. The density at the outer part of the 
atmosphere increased rapidly (over only 0.06 orbits) as a shock front propagated diagonally 
along magnetic field lines into the ray we are considering. The resulting changes in density 
structure as the shock passes through the ray results in the Fe Ka EW to vary drastically and 
show slight oscillatory behavior (-Fx/-^disk — 1)- We have thus clearly shown that changes 
in the accretion disk structure can cause rapid variability in the Fe Ka EW, despite the 
continuum remaining constant. 

The right-hand panel in Figure 1 illustrates the changes to the O VIII Lya line resulting 
from the evolution of the PBI. The line, of course, is weaker when the disk is more highly 
ionized at large i^x/^disk- The density enhancements that pass through the surface of the 
disk after about 1 orbit, decrease the local ionization parameter and increase the O VIII EW, 
often by > 10 eV, except for -Fx/-Fdisk =1 and 0.5. In these cases the increase in density makes 
it difficult for the weak continuum to ionize much oxygen, resulting in enhanced continuum 
absorption that suppresses that low energy reflection spectrum. In these cases, the O VIII 
EW will be (sometimes drastically) reduced. Unlike the Fe Ka line, the O VIII EW exhibits 
small amplitude rapid variability between orbits 1 and 2. This is because the soft X-ray 
lines are very sensitive to where the illuminated gas is most rapidly cooling. As a result, the 
strength of the O VIII line, as well as the N VII and C VI lines, are more sensitive to the 
density structure than the Fe Ka line (see also Paper I). 

Since the EW is a measure of the line strength relative to the continuum, and the illumi- 
nating continuum is held constant in these calculations, the variations in both the Fe Ka and 
O VIII EW are entirely due to intrinsic changes in the line flux in the reflection spectrum. 
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Therefore, measurements of the line flux would also see changes of similar magnitude due to 
this effect. 



3.2. Time Averaged Spectra 

The variations presented above occurred on timescales as short as a few one-hundredths 
of an orbit, or ~ 5-10 ks for the radius and black hole mass assumed for the simulation. 
This time is too short to accumulate a high quality spectrum with modern observatories 
such as XMM-Newton and Chandra except for the brightest Seyfert 1 galaxies. Typical 
observations need to be many times longer for proper spectral analysis, which will average 
over any variations in the emission lines due to the changing accretion disk structure. Thus, 
it is interesting to check whether the time-averaged spectra presented in § 3.1 retain any 
knowledge of the unstable disk structure. 

We follow very similar procedures as described in Paper I. For each value of i^x/^disk, 
we averaged together 156 reflection spectra between orbits 0.8 and 2.35 (in order to concen- 
trate on the region where the PBI is most important) and simulated a 40 ks XMM-Newton 
observation 1 with XSPEC v. 11. 3. Op (Arnaud 1996). The T = 2 power-law was added to 
each averaged reflection spectrum so that the reflection fraction R is unity. The simulated 
spectra were then fit with reflection models computed by the same code assuming a constant 
density of n H = 10 13 cm~ 3 . The results of the spectral fitting are shown in Table 1. The 
fits are generally acceptable in the 0.2-12 keV energy range, but are greatly improved in 
the 2-10 keV band. As a result of their sensitivity to the density structure, the soft X-ray 
emission features caused most of the difficulty in fitting the 0.2-12 keV simulated data (see 
also Paper I). The worst reduced x 2 obtained was for the i*x/-Fdisk=0.5 model over the 0.2- 
12 keV band. The residuals to this fit are shown in Figure 3, and are due to the slab model 
incorrectly accounting for the soft X-ray emission lines. However, the high energy continuum 
is well modeled by a uniform reflector. This exercise suggests that the density variations 
at the surface of the disk will not greatly affect the time-averaged reflection spectra above 
2 keV (but see § 4.2). Moreover, constant density reflection models will be useful to param- 
eterize observed spectra with an ionization parameter and reflection strength provided that 
the observing time is long compared to the fluctuation timescale in the disk. 



1 Strictly speaking, a 780 ks observation would be needed to observe a full orbit of this PBI simulation. 
However, since the normalization of the model is chosen so that Poisson noise never dominates, the actual 
length of the simulated observation does not change the results. 
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4. Discussion 

4.1. Variable Reflection Features 

In Paper I we showed that density inhomogeneities in the surface layers of radiation 
dominated accretion disks could potentially impact the reflection spectrum. That paper 
also suggested that the motion and evolution of any inhomogeneities in a real accretion disk 
might cause line variability that is independent of the illuminating radiation field. Here, we 
have confirmed that suggestion by explicitly showing variations in the Fe Ka and O VIII 
Lya EW on timescales as short as a few one-hundredths of an orbital time as a result of 
the evolution of the PBI, which is relevant for radiation dominated accretion disks. This 
mechanism provides a natural explanation for any short timescale Fe Ka variability that 
is uncorrelated with the continuum (e.g., Nandra et al. 1999; Wang, Wang & Zhou 2001). 
In particular, Fig. 1 indicates that the sudden change in Fe Ka EW observed in Mrk 841 
(Petrucci et al. 2002) could easily be accounted for by changes in the accretion disk structure. 
Furthermore, as mentioned in Paper I, the narrow, redshifted lines that have been recently 
inferred in the spectra of NGC 3516 (Turner et al. 2002) and Mrk 766 (Turner, Kraemer & 
Reeves 2004) could originate in density inhomogeneities at a particular radius in the disk 
atmosphere. If this is the case, then our results indicate the lines should be variable on 
timescales about 1% of the orbital period of the emitting patch. While the statistics are 
generally poor, there is evidence for variability in these narrow lines. Future data will be 
able to provide a better test of this hypothesis. 

The time-averaged reflection spectra from the PBI calculation showed no obvious indi- 
cation of the variations on smaller timescales. However, we cannot rule out rapid variations 
in disk structure that produce iron Ka changes drastic enough to be detected in observations 
lasting several orbits. In this paper, we considered only one black hole mass, accretion rate, 
position in the disk and magnetic field arrangement. For the parameters chosen, the Fe Ka 
region is readily fit with a uniform-slab reflection model, while the soft X-ray region in some 
cases is not. Clearly, there is wide parameter space that needs to be explored in order to 
determine the full impact the PBI and MRI will have on X-ray reflection spectra. 

The results in § 3 have important implications for upcoming Fe Ka reverberation 
mapping attempts with Constellation- X. If accretion disk structural changes can cause 
the Fe Ka line to vary independently of the continuum on timescales much smaller than 
£dyn ~ V8(R/Rs) 3 ^ 2 GM B}1 /c 3 (~ 1.2 x 10 5 s, for the photon bubble simulation shown here), 
then it will become an important source of noise for the reverberation experiments. This will 
be especially true for reverberation experiments probing small radii around low mass black 
holes (say, M B h ~ 10 6 ~ 7 M ). In that case, the disk structure may vary on timescales as 
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short as 100 seconds or less, comparable to the X-ray variability timescale (e.g., Vaughan, 
Fabian & Nandra 2003). Furthermore, the current transfer functions computed for reverber- 
ation mapping have all neglected the effects of ionization of the disk surface (e.g. Reynolds 
et al. 1999), which can result in an ant i- correlation between the line EW and continuum 
(Ballantyne & Ross 2002). Therefore, it is possible that the combined effects of a variable 
density structure (due to the PBI and perhaps other instabilities) and ionization state (due to 
changes in the X-ray illumination) will significantly complicate the results of any reverbera- 
tion experiment. What is needed to provide concrete predictions for reverberation mapping 
is time dependent ionized reflection calculations where the illuminating flux is related to 
the underlying magnetic stress (e.g., Armitage & Reynolds 2003) or outgoing magnetic flux 
(Miller & Stone 2000). This is a direction for future research. 

In the standard reverberation mapping scenario, the idea is to observe changes in the 
Fe Ka line, as X-rays from a flare sweep across the disk. The PBI will likely be ongoing at 
every radius in the inner region of the disk, but will be subject to local conditions, such as 
the strength and geometry of the magnetic field. Moreover, the linear growth rate of the 
PBI increases as one moves to smaller radii (Blaes & Socrates 2001). Therefore, the disk 
structure will be very position dependent as the X-rays flash overhead, so that the reflection 
spectrum could show variability due to disk inhomogeneities as a function of position, as 
well as time. PBI simulations spanning a larger range in radius would be necessary to test 
the possibility of reflection spectra variability due to positional changes. 

If different radii produce different variability properties, than the observed variations 
may be diluted if the X-rays illuminate a wide enough range of radii. However, there are 
often times when the AGN continuum is likely dominated by a single flare or a number of 
small flares which occur in close proximity to one another (cf. Merloni & Fabian 2001). In 
this case, which is the most relevant for reverberation mapping, then the one-zone models 
employed here may be useful. 

The signature of a changeable disk atmosphere is, of course, short timescale Fe Ka 
(and O VIII Lja, if observable) variability which is independent of the observed contin- 
uum. Therefore, it may be possible to observationally test if this process is important. As 
mentioned above, such cases of rapid changes in the Fe Ka line have been observed (e.g., 
Mrk 841), but these cannot be used as evidence that effects of disk instabilities have been 
observed (although it may provide a good explanation). There still remains too many un- 
knowns regarding the location and dynamics of the X-ray source providing the illuminating 
continuum. The long timescale disconnectedness between the Fe Ka line and continuum 
(e.g., Markowitz et al. 2003) suggests it is possible, if not likely, that at some level the 
disk sees a different irradiating continuum than what is observed. Numerical simulations 
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of coronal formation above accretion disks, still some years away, may provide the clues to 
understanding the X-ray source in AGN. 

4.2. Implications of the Time Averaged Results 

The reflection spectra that were averaged over the variability caused by the PBI were 
well described by a constant density model in the 2-10 keV band. The implication is that, 
over this energy range, the spectrum is dominated by emission from a single ionization 
parameter. Since the illuminating flux is known for these test cases, the best fit £ can 
be used to determine = 47ri 7 x/£, the density at which the most important spectral 
features are formed (see Table 1). The value of nu is always less than 10 13 cm -3 , the density 
actually used in the constant density models, indicating that the ionization level is higher 
than 'expected' given the illuminating flux. As can be seen in the density profiles plotted 
in Fig. 2, n H predominantly falls within the first 1-2 Thomson depths of the atmosphere. 
We conclude that reflection spectra above 2 keV are effectively determined by the ionization 
parameter at a Thomson depth of ~ 1. At lower energies, the carbon, nitrogen and oxygen 
cooling lines, which are very sensitive to the density structure, cause the spectra to exhibit 
features from different effective ionization parameters, and thus are more difficult to fit with 
a uniform density model. 

The fits with the constant density models suggest that changes in the calculated spectra 
on 0.01-0.1 orbit timescales are difficult to detect in the longer integrations. However, there 
may be regimes where rapid changes in disk structure produce clear signatures in the Fe 
Ka region in observations lasting several orbits. As mentioned in § 4.1, we considered only 
one black hole mass, accretion rate, position in the disk and magnetic field arrangement. 
Thus, there is a wide parameter space to be explored to determine the impact of accretion 
disk inhomogeneities on X-ray reflection spectra. Nevertheless, insofar as can be determined, 
these results (along with those from Paper I) indicate that constant density models provide 
a good indication of the average ionization state of the accretion disk atmosphere on long 
timescales. 



5. Conclusions 

Our conclusions from the second paper in our study on X-ray reflection from heteroge- 
neous accretion disks can be summarized as follows: 

1. Changes in the density structure of an accretion disk atmosphere from the photon bub- 
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ble instability can give rise to rapid variability of emission lines in the X-ray reflection 
spectrum. These variations would be observed to be independent of the illuminating 
continuum. 

2. This process gives a natural explanation for any short timescale (e.g., hours) line 
variability that seemed to be uncorrelated with the X-ray continuum (for example, 
Mrk 841). However, these observations cannot be used as proof that disk instabilities 
are causing the line variations. Nor can it offer an explanation for the uncorrelated 
behavior of the Fe Ka line on longer timescales (greater than a day), although we can 
not rule out such an effect in a different region of parameter space. 

3. The timescale of the variability was found to be as small as a few one-hundredths of 
an orbital time. Therefore, this effect, combined with any changes in ionization state, 
could be a serious source of noise for future reverberation studies. 

4. Averaging the individual reflection spectra over the 1.55 orbits most affected by the 
PBI yielded spectra corresponding to the more common long timescale X-ray AGN 
observations. These averaged spectra were well fit above 2 keV by constant density 
models. 

5. Photon bubble structures move across the surface of the disk in the radiation MHD 
calculations, passing in and out of the ray where we extracted density profiles. Varia- 
tions in the X-rays reflected from a region containing many structures may be less than 
those from a single line of sight. The separations, strengths and orientations of the 
structures are found to vary with the alignment and strength of the magnetic fields. 
The fields are likely to evolve over orbital timescales along with the turbulence inside 
the disk. Two and three-dimensional reflection calculations may be useful in studying 
the consequences for reverberation mapping. 
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Table 1. Results of fitting the time averaged photon bubble reflection spectra with constant density models. 







0.2-12 keV 






2.0-10 keV 






Fx/F disk 


r 


£ R 


xVd.o.f. 


r 


t R 


X 2 /d.o.f. 


n H 



0.5 2.012±0.002 126±0.3 1.0211^ 3832/2358 2.0171^ 741^ 0.8911^^ 1717/1595 6.5X10 1 

1 2.005i° ) ;g°J 249±10 0.9931°,;^ 2706/2358 2.008±0.010 219l|l 0.889^;}^ 1714/1595 4.4X10 1 

4 1-977+^ 6971}? 1.109±g;gfg 3814/2358 2.022±°;°}° 964+^° 0.856^;^ 1615/1595 4xl0 12 

8 1.997±g;ggf 1122t? 8 3 0.837+°,;^ 3432/2358 2.024±g;g5o 1738 -ii - 912 ±o:o49 1554/1595 4.45x10 



Note. - Averaged between 0.8 and 2.35 orbits. The normalization of the photon bubble models were chosen so that the 
'observed' 2-10 keV flux was 3-4xl0 -8 erg cm -2 s~b R is the reflection fraction defined as total=incident+i?x reflected, and 
d.o.f.=degrees of freedom. £ = 47rFx/nn is the ionization parameter and has units of erg cm s~h % is the density derived 
from the known Fx and fitted £. It has units of cm~ 3 . The error-bars are the 2a uncertainty for the parameter of interest. The 
simulated data were constructed using the XMM-Newton response matrix epn_sw20_sdY9 . rmf and assumed an exposure time of 
40 ks. 
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Fig. 1. — The evolution of the Fe Ka and VIII Lyct equivalent width (EW) as the photon 
bubble simulation runs over 2.35 orbits. The T = 2 power-law was added to the reflection 
spectra before the EW was calculated. The Fe Ka EW exhibits variability over timescales 
as small as a few one-hundredths of an orbit, and can change by factors up to 2.5. The large 
drop of the VIII EW when -F x /-F disk = 0.5 is due to increases in density causing significant 
absorption of the reflection spectrum. 
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Fig. 2. — X-ray reflection spectra (solid line) calculated from the outer 10 Thomson depths 
of a photon bubble simulation when Fx/F disk = 4. The short-dashed line denotes the T = 2 
power-law that was incident on the material. The insert shows the gas temperature (solid 
line) and number density (dashed line) for each model. The different plots show how the 
reflection spectra, temperature and density changed as the simulation evolved from Orbit 
1.6 to Orbit 1.69. 
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Fig. 3. — Simulated count spectra and residuals (in units of standard deviations) of the 0.2- 
12 keV constant density fit to the time averaged Fx/^disk =0.5 model. The fit parameters 
are given in Table 1. 



